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We establish estimates for the local and uniform moduli of continuity of the local time of 
multifractional Brownian motion, B H = {B H ^\t),t £ An analogue of Chung's law of the 
iterated logarithm is studied for B H and used to obtain the pointwise Holder exponent of the 
local time. A kind of local asymptotic self-similarity is proved to be satisfied by the local time 
of B H . 
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1. Introduction 

Multifractional Brownian motion (mBm), B H = (B H ^ (t),t £ is a Gaussian process 
which extends fractional Brownian motion (fBm) by allowing the Hurst function to vary 
with time. This provides a tool to model systems whose regularity evolves over time, such 
as internet traffic or images. Recently, Boufoussi et al. [6] have investigated, under mild 
regularity conditions on -ff(-), the existence of jointly continuous local times of mBm. 
Their objective was to explore the link between the pointwise regularity of the local 
time and the irregularity of the underling process. This effect of inverse regularity was 
observed by Berman [5] for the case of uniform regularity. 

The first aim of this paper is to establish estimates for the local and uniform moduli 
of continuity for the local time of mBm. Upper bounds for the moduli of continuity 
of local times have been obtained by Kono [16] for Gaussian processes with stationary 
increments and, more recently, Csorgo et al. [8] have proved upper bounds for the moduli 
of continuity of the maxima of local times for stationary Gaussian processes and Gaussian 
processes with stationary increments. However, contrary to classical results, the moduli of 
continuity obtained in this paper depend on the point at which the regularity is studied. 
This is natural, since mBm has a regularity evolving over time. The approach of the 
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present paper seems to be useful for extending the results of Kono [16] and Csorgo et al. 
[8] to Gaussian processes without stationary increments. 

Chung's form of the law of iterated logarithm for mBm is established in Section 4. 
This result is used to prove that the pointwisc Holder exponent of its local time with 
respect to time is equal to 1 — H(t), uniformly in x. 

The second main objective of the paper concerns the question of whether the sample 
path properties of mBm can be transferred to its local time. If H(-) is not constant, the 
mBm is no longer self-similar. However, it is proved by Levy-Vehel and Peltier ([17], 
Proposition 5) that if H is /?-H61der continuous and sup tgR + H(t) < (3, then a property 
called local asymptotic self- similarity remains, defined as follows: 

lim + lawi BiP+ p~ B ^\ n G R j = law{B H ^ («), iieM}. (1) 

where B H ^ is fBm with Hurst parameter H{t). It is the purpose of Section 5 to prove 
that the local time of mBm has a kind of local asymptotic self-similarity. Through this 
result, we obtain some local limit theorems corresponding to the mBm. 

We will use C, Ci,... to denote unspecified positive finite constants which may not 
necessarily be the same at each occurrence. 



2. Preliminaries 



In this section, we present some notation and collect facts about multifractional Brownian 
motion and its local times. 



2.1. Multifractional Brownian motion 

mBm was introduced independently by Levy-Vehel and Peltier [17] and Benassi et al. 
[3]. The definition due to Levy-Vehel and Peltier [17] is based on the moving average 
representation of fBm, where the constant Hurst parameter H is replaced by a function 
H(t) as follows: 

+ (i-w) ff(t) - 1/2 VF(dw)^, t>0, (2) 

where H (•) : [0, oo) — > [fi, v] C (0, 1) is a Holder continuous function and W is the stan- 
dard Brownian motion defined on (— oo,+oo). Benassi et al. [3] defined mBm by means 
of the harmonizable representation of fBm as follows: 

B H (Ht)=ljf^*W0, (3) 
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where W(£) is the Fourier transform of the series representation of white noise with 
respect to an orthonormal basis of L 2 (R) (we refer to Cohen [7] for the precise definition 
and for the fact that B H is real- valued) . 

Various properties of niBm have already been investigated in the literature, related, 
for instance, to its pointwise and uniform Holder regularity, as well as the local Hausdorff 
dimension of its sample paths. More precisely, it is known from Levy-Vehel and Peltier 
([17], Proposition 8) that with probability one, for each t > 0, the Holder exponent at 
to of mBm is H(ta). Recall that the pointwise Holder exponent of a stochastic process 
X at to is defined by 

f , X(t + p,w) - X(t ,uj) ] „ x 

ax(t ,u) = supja > 0, lim 1 ' ' -L L22_^ =0 |. (4) 

In addition, according to the same authors, the local Hausdorff dimension of the graph 
is 2 — min te [ a a H{t) almost surely for each interval [a, b] c R + . 

According to the previous results, the regularity of mBm depends on the regularity 
of H. Furthermore, by using Lemma 3.1 in Boufoussi et al. [6], we can easily prove 
that the irregularity of H implies the irregularity of mBm. More precisely, the points of 
discontinuity of H are also discontinuities of B (see Ayache [2], Proposition 1). 



2.2. Local times 

We end this section by briefly recalling some aspects of the theory of local times. For a 
comprehensive survey on local times of both random and non-random vector fields, we 
refer to Geman and Horowitz [12], Dozzi [10] and Xiao [21]. 

Let X = (X(t),t G R + ) be a real-valued separable random process with Borel sample 
functions. For any Borel set B C R + , the occupation measure of X on B is defined as 

Hb(A) = A{seB: X(s) G A}, for all A G B(R), 

where A is the one-dimensional Lebesgue measure on R + . If fis is absolutely continuous 
with respect to Lebesgue measure on R, we say that X has a local time on B and define 
its local time, L(B, •), to be the Radon-Nikodym derivative of [Ib- Here, x is the so-called 
space variable and B is the time variable. 

By standard monotone class arguments, one can deduce that the local times have a 
measurable modification that satisfies the following occupation density formula: for every 
Borel set B C R + and every measurable function / : R — > R + , 

J f(X(t))dt = J f(x)L(B,x)dx. 

Recently, Boufoussi et al. [6] have proved that if the Hurst function is Holder continuous 
with exponent (3 and if sup t>0 H (t) < (3, then the mBm has a jointly continuous local 
time, that is, the mapping (t,x) — > L{t,x) is continuous. In addition, this local time has 
the following Holder continuities. It satisfies, for any compact U C R, 
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(i) 

sup |L(t + ^,)-L(^)| <+oo 
xeu W 1 

where 7 < 1 — H(t) and \h\ < 77, r\ being a small random variable almost surely positive 
and finite; 

(ii) for any / C [0, T] with small length, 

\L(I,x)-L(I,y)\ 
sup _ < +00 a.s., (6) 

wherea<( 2supj 1 g(t) -§)Al. 

These results have been used to obtain the local and pointwisc Hausdorff dimension 
of the level sets of mBm. We refer to Boufoussi et al. [6] for definitions and results. 

3. Moduli of continuity of the local time 

Throughout this section, the Hurst function H — > [fj,,v] C (0,1) is assumed to be 
measurable. The notation B = (B(t),t > 0) means that either representation, (2) or (3), 
can be chosen. Moreover, we say that H satisfies the condition (Tip) if 

H is /3-H61dcr continuous with sup-ff (s) < (3. 

s>0 

We give the following improvement of Theorem 3.1 of Boufoussi et al. [6]. Note that 
the first part reproduces this theorem without assuming {Tip). 

Theorem 3.1. Consider a measurable function iJ(-):R_|_ — -> \p,,v] C (0,1). The mBm 
with Hurst function H(-) given by (2) admits on any interval [a,b] C [0, 00) a square- 
integrable local time. Moreover, for both representations (2) and (3), if (Hp) holds, the 
existence of square-integrable local times implies that H(t) < 1 for almost all t. 

Proof. Let us write, for simplicity, 

B H ^(t)= f K H(t) (t,u)dW(u), 
J —00 

where 

For any t > s, taking the variance of B H ^'(t) — B H ( s '(s), we obtain 

\&y(B H ^ (t) - B H ^ (a)) > Vax(B H ® (t) - B H ^ (s)/W(u),u < s) 

= V&r(B H ^(t)/W(u),u<s), (7) 
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where the last equality follows from the fact that B H ^ (s) is measurable with respect to 
a(W(u), u < s). Moreover, we can write 

B H &(t)= [ K H{t) (t,u)dW(u) = f K H(t) (t,u)dW(u) + f K H(t) (t,u)dW(u). 



Hence, by using the measurability of KhU) (t, u) dW(u) with respect to a(W(u)/u < 
s), we obtain 

V&T(B H W(t)/W(u),u< s) =Vax(J^ K H(t) (t,u) dW(u)/W(u),u < s^j 

= Var(7 K m (t,u)dW(u)\ (8) 

where, to obtain the last equality we have used the fact that J K[j(t){t, u) dW(u) is 
independent of <j(W(u),u < s) (by the independence of the increments of the Brownian 
motion). Combining (7) and (8) and denoting C = sup u6 [ M „j r(l/2 + u), we obtain 

Var(B ff «(t)-B"«( S )) > ± j\t - rf H ^ dr 

> 1 (t-s) 2Hit) 



Therefore, 



/ / (E[B H ^(t)-B H ^( S )} 2 )~ 1/2 dsdt 

J\a,b] J\a,b] 



l[a,b] J[a,b] 

<V^C [ [ \t~s\- sup r>" H ^dsdt. (9) 

J[a,b] J[a,b] 

The last integral is finite because sup r>0 H{r) < 1. By Theorem 22.1 in Gcman and 
Horowitz [12], B H has local time L([a,b], •) e L 2 (R). 

We now prove the second point. According to Boufoussi et al. [6], we have, for any 
interval [a, b] C [0,oo) with small length, 

E(B(t)-B(s)f<C^\t-s\ 2H ^, for s,t€[a,b}. (10) 

On the other hand, since the local time exists on the interval [a, b], then according to 
Geman and Horowitz ([12], Theorem 22.1 expression (22.3)), the following integral is 
finite 
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Then, (10) implies that 



6 r b 

ds dt < 



la J a (t-s)HW 

Consequently, H(t) < 1 for almost all t £ [a, b]. Since R+ is a countable union of small 
intervals, the result is proved. □ 

In what follows, we are interested in the local and pointwisc oscillations (at each t) of 
the local time of mBm. 

Theorem 3.2. Let {B H{ $ (t),t > 0} be the mBm given in (2) with Hurst function H(-) 
satisfying the assumption (Hp). Then, for every t £ M. + and any x £M., there exist positive 
and finite constants C\ and C2 such that 

v L(t + 5,B H V(t))-L(t,B H M(t)) <r 

1 T-! up gi-g('>(io S b 8 (g-i))g(') ~ Cl a ' s -' 

5^0 S 1 - H W(loglog(S- 1 )) H< - t > 

Proof. Let t > be a fixed point, suppose < h < 1 and define H t: t+h = su P s e[t,t+W H( s )- 
According to the Fourier analytic approach of Berman [4] (see also Davies [9], expres- 
sion (27)), we have 

E[L(t + h, B H ^ (t)) - L(t, B H ^ (t))] m 

Let B(s) = B H ^ (s) ~ B H ^ (t), s> 0, and denote by R(s\, . . . , s m ) the covariance matrix 
of (B(s 1 ), . . . ,B(s m )) for distinct si, . . . ,s m . Let U = (ui, . . . , u m ) £ R m and let U T de- 
note the transpose of U. According to (16) below, we have deti?(si, . . . , s m ) > 0. Hence, 
the change of variable V = R}I 2 U implies that 



f E(e i ^r= 1 "^ ) )d Ul ---d^ 



(27r)™/ 2 



(deti?( Sl ,...,s m ))V2 
Therefore, 

E[L(t + h, B H ^ (t)) - L(t, B H ^ (t))] m 

>/* J [t , t+h]m (deti?( Sl ,..., Sm ))V2 dsi ' ' ' ds '" 



(27r)™/2 J t<Sl< ... <Sm<t+h (detR( Sl , . . . , Sm ))V 



jTjdsi--- ds m . (13) 
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deti?(si, ...,s m ) 

= Var(B( Sl )) Var(5(« 3 )/B(«i)) ■ ■ • Var(B(s m )/5( Sl ), . . . , S(s TO _i)). (14) 

Using arguments similar to those used in the proof of Theorem 3.1, wc obtain, for any 
r,s£ [t,t + h] such that r < s, 

Var(B(s) /B{u), ueA,u<r)> Var(£?(s) - B(r)/W(u), u < r) 

> Var(S ff(s) (s)/W(u), u < r) 

where C = sup^r^i T(.t + 1/2). Combining (14) and (15), we obtain 
deti?( Sl , . . . , s m ) > H( Sj - 8j -x) aH M 

^j^U^-^r^, da) 

where s = and (sj — Sj-i) H ^ > (sj — Sj-\) Ht < t + h , since (sj — Sj-i) < 1. 
According to (13), we have 

E[L(t + ft, B H{t) (t)) - L(t, (t))] m 

"V T J m '7t< Sl <...< Sm<t+ / l /J («j r-B i _i)*M+* 

Now, by an elementary calculation (cf. Ehm [11]), for all m > 1, ft. > and bj < 1, 



Jt< Sl <—<s m <t+h - =1 1 V 1 T"l 2^j=l°3J 



Therefore, 

£[L(t + ft, (t)) - L(f , (t))] 
<m\r—\ h 



T(l + m(l-H t , t+h ))- 
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According to Stirling's formula, we have (m!/r(l+m(l-F M+h ))) < M m m\ H ^+ h ,m> 2, 
for a suitable finite number M . Therefore, 



E 



L(t + h, B H ^ (t)) - Lit, B H ^ (t)) 



< C m m\ Ht t + h . 



(17) 



We shall now prove that for any K > 0, there exists a positive and finite constant A > 0, 
depending on t, such that for sufficiently small u, 



P \L(t + h, (t)) - L(t, B*<«> (t)) > ugtt+h j < exp(-K/ U ). 



(18) 



First consider it of the form u = 1/m. Combining Chebyshev's inequality and (17), we 
obtain 



L(t + h,B H M(t))-L{t,B H M(t)) 



<E 



fol — H t; t+h 

L(t + h, B H ® ( t )) - L(t, B H ^ (t)) 



> Am Ht t+h 



< 



Ah^ 1 ^ Ht - t+h ^m Ht - t + h 



A m \m 



(ml) 



Again using Stirling's formula, the last expression is at most ^ T (2Trra) Ht - t + h / 2 c H t,t+h m . 
This can be written as 



m[log(C/A) - H t , t+h ) + ^±[l g(m) + Iog(27r)] 



(19) 



Choose A> C and mo large such that for any m > mo, to dominate (19) by e~ 2Km . 
Moreover, for u sufficiently small, there exists m > mo such that u m +i < u < u rn and 
since m > 1, > |. This proves (18). 

In addition, if we take u(h) = l/loglog(l//i) and first consider h m of the form 2~ m , 
then (18) implies 

P(L(t + h m , B H ^ (t)) - L(t, B H ^ (t)) > Ah l - Ht > t+h ™ {\og\og(l/h m )) Ht ' t+h ™ ) < m- 2 

for large m. Consequently, by using the Boref-Cantelli lemma and monotonicity argu- 
ments, we obtain 



L(t + h, B H ^ (t)) - L(t, B H ^ (*)) 



<A(loglog(l//i)) 



Ht,t + h 



In addition, the Holder regularity of H implies that 

loglog(l/fe)^ = Kh e logloglog 1/h 
loglog(l//i) ff(t) 



a.s. 
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Therefore, 

(Ioglog(l/fc)) H *^ < K(loglog(l/h)) H(t) 

for sufficiently small h. Moreover, using the same arguments (see Remark 3.7 in Boufoussi 
et al. [6]), we obtain 

h i-H t , t+h <jf ft i-ff(t). (20) 

This completes the proof of (11). Since (12) is proved in the same manner, the proof is 
omitted here. □ 

Remark 3. 3. The result of the previous theorem, with sup sg u t+h ^ H(s) instead of H (t) , 
may be proven even when (Tip) is not satisfied. 

We shall also prove the following uniform result. 

Theorem 3.4. Let {B H ^(t),t € [0, 1]} be the mBm in (2) with arbitrary Hurst function 
and define H i = sup sg [ H(s). Then, for all x G R, there exists C3 > such that 

limSU Pl J L M ,X n~^n ,X)l <C 3 0.8. (21) 

t,ae[0,l] 

Proof. By using the same arguments as above, without using the Holder continuity of 
H , we show that with probability one, there exists mo(w) such that for all m > mo(w) 
and fc=l,...,2 m , we have 

L(k2- m ,x) - L({k - l)2~ rn ,x) < A2- {1 - Ha ^ m (m\og{2)) Ha ' 1 a.s. 

By then proceeding essentially as in Kono [16], we prove (21). □ 

Remark 3.5. The results in this section are based on the following property of the 
mBm: 

Var(S ff {t+h) (t + h) /B H{S) (s):0<s<t)> Ch 2H ^+ h , (22) 

for sufficiently small h. If H is constant, this property is called one-sided strong local 
nondeterminism; we refer to Monrad and Rootzcn [20], expression (2.2), and the refer- 
ences therein for definition and applications. Property (22) is satisfied by many other 
processes with multifractal behavior, such as the Riemann-Liouville mBm {X(t),t> 0} 
introduced in Lim [18] by replacing, in the definition of the Riemann-Liouville fBm, the 
constant Hurst parameter H by a regular function as follows: 
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Thus our results extend those of Kono [16] to Gaussian processes without the stationary 
increments property. 

Remark 3.6. A closer examination of the proof of Theorem 2.1 in Csorgo et al. [8] 
shows that the property of stationarity of increments is only used to prove two main 
ingredients presented in their Lemmas 3.1 and 3.3 and which can, in our case, be easily 
replaced by the following properties: 

1. Let A n be the covariance matrix of a Gaussian vector {Q; 1 < i < n}. It is then 
known that the conditional variance can be written as 



Vax(Q/Q,l^i,l<l<n) = 




where A„ is the submatrix of A n obtained by deleting the ith row and column. 
2. If a stochastic process has the property (22), with H constant, then, according to 
(14), its covariance matrix satisfies 

m 

det R(t u . . .,t m ) > C m Y[(tj-tj-i) 2H , for t <*!<••• < t m <t+h and t Q = 0. 

3 = 1 

The arguments of Csorgo et al. ([8], Theorem 2.1) can then be adapted to prove that the 
class of processes satisfying (22) satisfies 

v L(t + a,T,x) —L(t,x) 

hmsup sup — j— 77 ■ — < C < oo a.s., (23) 

T-+oo o<t<b T a T (\og(b T /aT) + loglog(aT + 1/ar)) 

where ax and &t are non- negative functions in T > such that 1 "^ bT — > +oo as T — > +oo. 



4. Chung's law of the iterated logarithm 

The following theorem establishes an analogue of Chung's law of the iterated logarithm 
for mBm. The result for fBm has been obtained by Monrad and Rootzen [20]. 

Theorem 4.1. Let B be a multifractional Brownian motion and assume that H satisfies 
condition (Tip). Then the following Chung-type law of iterated logarithm holds: 



\B(s)-B(t ) 
^°\e[wU] (*/log|bg(*)|)*(M 



limi ? f SU P ^'■^L/^ffL =C H (t ) a.s., (24) 



where the constant Cjj{t ) * s the one appearing in Chung's law for the fBm of Hurst 
parameter H(to). 
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Proof. We first introduce the process 0(t) = B(t) - B Hi - to \t),t > 0}, where B H ^ 
denotes a ffim with Hurst parameter H(to). According to Theorem 3.3 of Monrad and 
Rootzen [20], the fBm with Hurst exponent H(to) satisfies (24). Furthermore, 

B{t) - B(t ) = (B(t) - B H ^ tn \t)) + {B H{t °\t) - B H ^(t )), 

= B(t) + {B H ^(t) - B Hit "\t a )). 

(24) will then be proven if we show that 

\B(t)\ 

hm sup — — — r^KWn \ = a - s - 

On the other hand, Lemma 3.1 in Boufoussi et al. [6] implies that there exists a positive 
constant K such that 

sup E(B(t)) 2 < KS 213 . 

t£[to,t Q +6] 

Hence, according to Theorem 2.1 in Adler ([1], page 43) and a symmetry argument, we 
have 

P[ sup \B(t)\ >u) < 2P[ sup B(t)>u 

( (u-E(su Pte[to , to+s] B(t))f 
<4exp(^ ^_ ). 

For the sake of simplicity, let A = sup tg j to ta +s] By (25), we obtain 

E(A)< ( p( sup \B(t)\>x]dx 
Jo \te[t a ,t +s] J 



\te[t ,t +5] 

, + °° ( \x-EA} 2 , 
CX P1" KS 2 ) dx 



e~y / 2 dy 

-J2/KEK/5P 



It follows that 



(u - EA) 2 > \u 2 - (EA) 2 > \u 2 - 16Kn6 213 
Consequently, (25) becomes 

J2/3")' 



sup \B(t)\>u)<Ce X p[—^ s ). (26) 



t<E[t ,t +8] 
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Since H(t ) < (3, there exists < £ < /3 - #(t ). Consider 5„ = n V(2(e+H(to)-/3)) and 
u n = Sn ^- Therefore, according to (26), we have 

y^P\ sup \B(t)\ > u n ) < ^exp( -— n ] < oo. 

t^i \te[t ,t +8 n ] J ^ V A / 

It follows from the Borel-Cantelli lemma that there exists n = n Q (oj) such that for all 
n > no, sup se j to t Q +5 n i |-B(s)| < Sn^ ^^ almost surely. Furthermore, for <5„ + i < S < 5 n , 
we have 



sup \B{s)\ < sup |B(s)| 
se[t ,t +s] se[t ,t +s„] 

— u n — 

< 2 e 8 H ^+t a.s, 

2(/3-ff(t )-£) 



where 6 = *J$+* . Keace, 



Urn sup — iM— — < lim tf«(log | log(S)\) Hito) = 



a.s. 



s ^°te[t ,t +5] («5/log|log(5)|) J 
This completes the proof of the theorem. □ 

Remark 1^.2. Observe that by the ideas used to prove the previous theorem, many laws 
of iterated logarithm (LIL) proved for fBm can now be obtained for mBm and with the 
same constants. For example, we have the LIL (cf. Li and Shao [19], equation (7.5) for 
the fBm) 

limsup sup \f ^ 7w»! u/ , = V2V H(to) (B) a.s., (27) 

5-0 s6 [t ,t +<5] ^^HlogllogWI) 1 / 2 

where 



H(t )T(2H{to))sm(irH(t )) 
and 

_ (f°J(l - s) mto)-l/2 _ ( _ s) H(to)-V 2]2ds) l/2 + l/ {2 H(t )) 

VH{t0)[B ' T(H(t ) + l/2) ' 



Other laws of iterated logarithm and uniform moduli of continuity have been obtained 
for B by Benassi et al. [3], Theorem 1.7, via wavelet techniques. 



Sample path properties of the local time of multifractional Brownian motion 



861 



Since the local time is a two-parameter process, we make precise the meaning of the 
pointwisc Hlder exponent in time, uniformly in space. The pointwise Hlder exponent of 
sup x&R L(-,x) is defined by 

M f ^ n i- L(t + S,x)- L(t ,x) \ 
OLL\t)= sup< a > 0,hmsupsup = >. (28) 

Since the fraction above is positive, the limsup may be replaced by lim or liminf. 
We are now able to prove the following. 



Proposition 4.3. Let B be a multifractional Brownian motion and assume that H satis- 
fies the condition (Tip). The pointwise Holder exponent ah o/sup l£ |,L(-,a;) at t satisfies 

a L (t) = l-H(t) a.s. 

Proof. The lower bound was already given in Boufoussi et al. [6], Corollary 3.6. The 
upper bound is a consequence of Chung's law of the mBm. Indeed, since the local time 
vanishes outside the range of B, we obtain 

5=1 L([t ,t Q + S},x)dx 
< sup L([t , t + 5], x) sup \B{s)-B(t)\ 

x£R s,te[t ,t +S] 

<2sup L([t ,t + 5],x) sup \B(s) - B(t )\. (29) 

xGR se[t ,t +S] 

Combining (24) and (29), we obtain that there exists a positive constant C such that 

,. L(t + 6,x) - L(t ,x) 

hm sup sup . — - — — ; > C a.s., 

5^0 ^eK^^^Hioglog^^ 1 ))^* ) 

which, together with the definition of the pointwise Holder exponent, proves the result. □ 



5. Asymptotic results 

It is well known that techniques for proving limit theorems related to self-similar processes 
use the self-similarity of their local times. It is natural to expect the same when dealing 
with locally asymptotically self-similar processes (LASS for short). Thus it will be of 
some interest to know if the local times satisfy a kind of LASS property. 



5.1. LASS for local times 

The answer to the preceding question is affirmative in the case of the mBm and the result 
is given by the following theorem. 
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Theorem 5.1. Let B be a multifractional Brownian motion and assume that H satisfies 
the condition (Ti-p). Then, for any fixed to, the local time ofmBm is locally asymptotically 
self-similar with parameter 1 — H(to), in the sense that for every i£R, the processes 
{Y p (t,x),te [0, 1]} P>0 , defined by 

v( , , _ L(t + pt, p H ^x + B(t )) - L(t , p H( - to) x + B(t )) 
r P^> x >- pi-ff(to) 

converge in law to the local time, {i(t,x),t G [0,1]}, of the fBm B Ht ^ to ^ with Hurst pa- 
rameter H(ta), that is, 

lim law{YJt, x), t G [0, 1]} = law{£(t, x), t G [0, 1]}, (30) 

p->0 

where the convergence is in the space of continuous functions endowed with the norm of 
uniform convergence. 

The motivation for considering the processes Y p follows from the occupation density 
formula and the LASS property (1) of mBm. Indeed, according to the occupation density 
formula, the local time has the following representation: 

L(t, x) = ^- [ + c~ iux ( [ c iuB[s) As\ du. 



2 71 

Consequently, 

1 f ■...r„H(t n )„, T1 u.\-\ f to+pt 



Y p (t,x) 



p pto-rpt 

Jr Jt 



Using the change of variables r = s — p ^ s - and p H ^y = v, the right-hand side of the previous 
expression becomes 

which is the local time of the Gaussian process {B' } {r) = B(t ° + p p H r { \~ } B(to) , r G [0, 1]}. We 
need the following lemma for the proof of finite-dimensional convergence. 

Lemma 5.1. Let {X(t), t G [0, 1]} be a stochastic process in the Skorohod space D([0, 1]) 
and define, for a fixed interval L = [a, b] C ffi, the map 

M*)(t)= I l{x (s)6 /}d S . 
Jo 

If {X n (t), t G [0, l]} ra >i is a family of processes which converges in law in D([0, 1]) to X , 
then 4>i{X n ) converges in law to 4>i(X). 
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Proof. The lemma is a consequence of the continuity of the map X(-) — ► (j>i(X(-)) in J\ 
topology on D([0, 1]) at almost all sample points of the process X, which is proved on 
page 11 in Kesten and Spitzer [15]. □ 

Proof of Theorem 5.1. To prove the convergence in law, we proceed in two steps. First, 
we prove the tightness of the family {Y p (t,x),t £ [0, l]} P >o in the space of continuous 
functions. By using (17) and (20), for sufficiently small p, we obtain 

E\Y p (t,x)-Y p (s,x)\ m 

_ E[L(t + Rt, p H ^x + B(t )) - L(t + ps, p H ^x + B(t Q ))] m 



p {l-H{t ))m 



<C m \t- s\ {1 - H ^ to))m . 



We can take m > 1 _^ to - j to prove the tightness. 

Now, we prove the convergence of the finite-dimensional distributions of Y p , as p tends 
to 0, to those of the local time t of the fBm B H ^ with Hurst parameter H(t ). We need 
to show that for any d> 1, <n, . . . , G M and t\ , . . . , G [0,1], the following convergence 
holds: 

d d 

^ajY p {tj,x)=^^ajl{tj,x) as p^O. 

3=1 3=1 

We will show the convergence of the corresponding characteristic function. More precisely, 
we will prove that 



Eexp 



i\J2 ajYpit^x) 

. 3 = 1 



Eexp 



. 3 = 1 



as p — ► 0. 



We introduce the following notation: 



0«,x (*)(*) 
1 1 



1 1 



Eexp 



Eexp 



i\^a 3 Y p {tj,x) 

L 3 = 1 



Eexp 



L 3 = 1 



L 3=1 



Eexp 



3=1 



and 



Eexp 



iAVa^(B^))(^) 



Eexp 



iA^~ ]a j£(tj,x) 
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Therefore, 



Eexp 



I 3=1 



Eexp 



L 3=1 



(31) 



On the other hand, Y p is the local time of B p and by using the mean value theorem and 
the occupation density formula, we obtain 



if' 



<C max E\Y p (t i ,x)-<t> e 0[B''){tj)\ 

l<j<d 



= C max E 

\<j<d 



Y p{tj-,v) dy-Y p (tj,x) 



(32) 



Moreover, since the stochastic process {Y p (t,y),y S M} is almost surely continuous in y 
for every t, according to the dominated convergence theorem, (32) converges to zero as 
£ tends to zero independently of p. 

We now deal with i|' p . Since the family of processes {B p (t),t£ [0, l]} P >o converges 
in distribution to the fBm {B H ^ to \t),t € [0, 1]} with Hurst parameter H(to), the second 
term converges to zero as p tends to by Lemma 5.1. 

The last term in (31) is treated in a similar way as the first and the proof of the 
finite-dimensional convergence is complete. □ 



5.2. Limit theorems 

The following result is an immediate consequence of the LASS property of mBm. 

Proposition 5.2. Let B be a multifractional Brownian motion, assume that H satisfies 
the condition {Tig) and denote by £(t,x) the local time of the fBm with Hurst parameter 
H(to). Then, for every f £ L 1 (K) which is locally Riemann integrable, with compact 
support and such that J R f(x) dx ^ 0, the following convergence in law holds: 

Um Hm [ Xt f( B{pS + t °l- Bito) )ds= [ f(x)dx.£(t,0). 

If H is constant, then B is fBm. Being self-similar with stationary increments, the term 
in the left-hand side of the previous formula is identical in law to jt^-ft J q M f{B(s)) ds. 
We retrieve the result of Kasahara and Kosugi [14] for the fBm. 

Proof of Proposition 5.2. The proof relies on the LASS property and the self- 
similarity of fBm. We sketch it for the sake of completeness. 

Combining the fact that / is locally Riemann integrable, the LASS property of mBm 
and Theorem VI. 4. 2 in Gihman and Skorohod [13], we obtain 

f / ( B( " a+ ;°, ) j B "° ) ) d. ^ J" / cfl-*i (.» ... » ^ o + . 
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Moreover, using the occupation density formula and the self-similarity of fBm, we obtain 

tj 57*77 f(B H{to) (s))ds^ [ f(x)£(t,\- H ^x)dx, 
* ' Ja Jr 

where = denotes equality in distribution. The result of the theorem then follows from 
the continuity of £{t, x,u) with respect to the space variable. □ 

However, the presence of a double limit may be not convenient and the result may be 
better if A and p are dependent. We prove the following local limit theorem. 



Theorem 5.3. Let B be a multifractional Brownian motion and assume that H satisfies 
the condition (Tip). Denote by £(t,x) the local time of the fBm with Hurst parameter 
H(to) and consider f G i 1 (R) such that 

[ \f(x)\ \x\^ dx < oo for some <£< . . — - and [ f[x) dx ^ 0.(33) 

Jm 2sup t>0 ii(i) 2 J R 



The following convergence in law then holds 

1 f pt+t0 jB(s)-B(t )-p H ^y\ J ( d) 



f 



0(P) 



ds 



f(x)dx-£(t,y), asp->0+, 



where 9(-) and </>(•) satisfy = p^M and _HpL_ =o(1) . 

Proof. Using the occupation density formula, we obtain 
1 f ( B(s)-B(t )-p H My ] 



*P(p) Jt 



0(P) 

L(pt + t o ,0(p)x + p H ^y + B(t )) - L(t , 6(p)x + p H My + B(t )) & 



■A.-H{t ) 



L(pt + t ,p H ^y + B(t )) - L(t ,p H ^y + B(t )) 



x / f(x)dx 



L(I ,d(p)x + p H ^y + B(t )) - L(I ,p H My + B(t Q )) ^ 



l-ff(to) 



(34) 
(35) 



where, in the last expression, we let Iq = [to, to + pt] for simplicity. According to The- 
orem 5.1, the expression (34) converges in distribution to £(t,y) x J R f(x)dx. It now 
suffices to prove that (35) converges to in some strong sense. We have 



E 



tt , L(I o ,0(P> + P Hito) V + B(t )) - L(I , p H ^y + B(t )) 
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< / \m\ 



L(I , 9(p)x + p H ^y + B(t )) - L(I 0:P H ^y + B(t )) 



P 



l-ff(to) 



dx.(36) 



L 2 (Q) 



Using essentially the same arguments as those used in proving expression (22) and Re- 
mark 3.7 in Boufoussi et al. [6], but for the process X(t) = B(t) — B(to) instead of the 
mBm, we obtain 



L(I , 6{p)x + p H ^y + B(to)) - L(I ,p H ^y + B(t )) 



P 



l--ff(to) 



L 2 (fi) 



for sufhcently small p and all < £ < 2sup 1 H ^ — \ . Hence, (36) is dominated by 



(37) 



This last integral is finite by assumption (33) and then (37) tends to zero as p tends to 
zero. This completes the proof. □ 
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